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Phase transition into the phase with Bose-Einstein (BE) condensate in the two-band Bose-Hubbard model 
with the particle hopping in the excited band only is investigated. Instability connected with such a transition 
(which appears at excitation energies 5 < \t'„\, where |/(',| is the particle hopping parameter) is considered. 
The re-entrant behaviour of spinodales is revealed in the hard-core boson limit in the region of positive values 
of chemical potential. It is found that the order of the phase transition undergoes a change in this case and 
becomes the first one; the re-entrant transition into the normal phase does not take place in reality. First order 
phase transitions also exist at negative values of 5 (under the condition S > 5crit ~ -0.l2|tQ|). At /.t < 
the phase transition mostly remains to be of the second order. The behaviour of the BE-condensate order 
parameter is analyzed, the (©,/i) and (|toliM) phase diagrams are built and localizations of tricritical points 
are established. The conditions are found at which the separation on the normal phase and the phase with 
the BE condensate takes place. 
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1. Introduction 

During the recent years Bose-Hubbard model (BHM) is proved to be a valuable tool in the 
theory of systems of strongly correlated particles. The model achieves a wide recognition due to a 
successful description of thermodynamics and dynamics of ultracold Bose atoms in optical lattices 
where a phase transition to the phase with the Bose-Einstein (BE) condensate (so-called Mott 
insulator (MI) - superfluid state (SF) transition) occurs at very low temperatures. Experimental 
evidence of BE condensation in optical lattices was found for the first time in works Q while 
theoretical predictions of such an effect were given earlier Q . Starting from the 90-ies of the past 
century, a series of papers were devoted to the theory of this phenomenon. Among the first key 
articles on the subject one should mention the work 0| where BHM was studied in the mean 
field approximation. The calculated therein phase diagrams demonstrate that in the simplest case 
(i.e., hopping of Bose particles in the presence of a single-site Hubbard repulsion) the MI-SF 
transition is of the second order. Moreover, it is supposed that particles reside in the ground 
state of local potential wells in the lattice. Forthcoming theoretical investigations in this field were 
performed with the use of various techniques, e.g., the random phase approximation (RPA) in the 
Green function method 0, 6|, a strong-coupling perturbation theory [3, Q, the dynamical mean 
field theory (Bose-DMFT) 9,[l3| as well as quantum Monte-Carlo calculations [HI, and other 
numerical methods. 

The Bose-Hubbard model is also intensively used for a theoretical description of a wide range 
of phenomena: quantum delocalization of hydrogen atoms adsorbed on the surface of transition 
metals l^, M 1 1 quantum diffusion of light particles on the surface or in the bulk (isl . , thermo- 



dynamics of the impurity ion intercalation into semiconductors |17l Il8l | . 

In the last mentioned applications, there is usually a restriction on the position occupation 
number (n^ ^1), which corresponds to the limit of an infinite Hubbard repulsion for the consid- 
ered model. Such a model of hard-core ions (where particles are described by the Pauli statistics) 
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is also known as the fundamental one for the investigation of a wide range of problems, e.g., super- 
conductivity due to a local electron pairing [19] or ionic hopping in ionic (superionic) conductors 

The study of a quantum delocalization or diffusion reveals an important role of excited vibra- 
tional states of particles (ions) in localized (interstitial) positions with a much higher probability of 
ion hopping between them [15l |22|, |23j . A similar issue of a possible BE condensation in the excited 
bands in optical lattices is also considered but the condition of their sufficient occupation due to 
the optical pumping (see, e.g. [24]) is imposed. An orbital degeneration of the excited p-state is 
accompanied by anisotropy of hopping parameters and causes the appearance of variously polar- 
ized bands in the one-particle spectrum. Such bands correspond by convention to different sorts 
(so-called "flavours") of bosons and their number correlates with the lattice dimensionality. In the 
framework of the necessary generalization of the Bose-Hubbard model, a possibility of the MI-SF 
transition to the phase with BE condensate in the pumping-induced quasi-equilibrium long-living 
state of the system has been established [25]. 

In the equilibrium case, the issue of BE condensation involving the excited states in the frame- 
work of ordinary Bose-Hubbard model was not considered in practice. The exception is the system 
of spin-1 bosons [26„ .27] where a hyperfine splitting gives rise to multiplets of local states resulting 
in closely-spaced excited levels. As demonstrated in (28l. [29|. the MI-SF phase transition could be 
of the first order when a single-site spin interaction is of the "antiferromagnetic" type. A similar 
change of the phase transition order also takes place for multicomponent Bose systems in the 
optical lattices (30| . 

In the present work we consider an equilibrium thermodynamics of the Bose-Hubbard model 
taking into account only one nondegenerated excited state on the lattice site besides the ground 
one. On the one hand, such a model corresponds to ID or strongly anisotropic (quasi-lD) optical 
lattice, and on the other hand, it is close to a situation that is characteristic of a system of light 
particles adsorbed on the metal surface. For example, the excited states of hydrogen atoms on the 
Ni(lll) surface are sufficiently distant so only the lowest one could be taken into account. We 
shall investigate a condition of instability of a normal state of the Bose system with respect to BE 
condensation considering a criterion of divergence of the susceptibility (x ~ {{ci\Cp)) ui\q=o.uj=o) 
characterizing the system response with respect to the field related to a spontaneous creation or 
annihilation of particles. We shall also study the behaviour of the order parameter (cq) {{cq}) as 
well as the grand canonical potential in the region of the MI-SF transition and shall build relevant 
phase diagrams. Special attention will be paid to a change of the phase transition order and 
localization of tricritical points at different values of excitation energy, particle hopping parameter 
and temperature. 

We shall limit ourselves to the hard-core boson (HCB) limit where a limitation on occupation 
numbers is present: no more than one particle per site regardless of the state (excited or ground) 
occupied by it. Thus, the single-site problem is a three-level one (contrary to the two-level ordi- 
nary HCB case) . For this reason, it is convenient to use the formalism of Hubbard operators (3l| 
(standard basis operators (32|). 



2. Two-state Bose-Hubbard model in RPA: normal phase 

The Bose-Hubbard model is used for description of the system of Bose particles which are 
located in a periodic field and can reside in lattice sites. Taking into account only the ground 
and the first excited vibrational levels in the potential well on the site, one can express the model 
Hamiltonian as: 

H = is-^,)J2 Kb^ + is' - A^) E + Y E "'(^' - 1) + Y E <« - 1) + E 

i i i i i 

+ E ^^^tb, + E ^A^^ + E ^^^^ + ^tb,), (2.1) 

ij ij ij 
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where bi and bf {ci and cf) are Bose operators of annihilation and creation of particles in the 
ground (excited) state, e and e' are respective energies of state and /i is the chemical potential 
of particles. Such a Hamiltonian includes the single-site Hubbard repulsions with energies Ub^ Uc 
and Ubc as well as the particle hopping between ground (t^), excited (t"^) and different (i'"^) states. 
Hereinafter we assume Ub = Uc = Ubc for simplicity. 

Let us define a single-site basis |n,^,n^) = |i;n^,nj^) (which is formed by particle occupation 
numbers in the ground and in the excited states, i.e., eigenvalues of operators = 6^6^ and 

= cfci) as well as introduce Hubbard operators (standard basis operators) 

A J ' ' ' =\i;n,m){i]n ,m\. (2.2) 
Annihilation and creation Bose operators may be written as 

bi = 2^2^Vn + lX.' • • , 6+ = ^ ^ + 1 ' ' ; 

71 in n rn 

Q = 2^2^ Vm-MX^ ■ ' , c+ = 2^ 2^ Vm + IX. ' . (2.3) 

n m n rn 

Corresponding occupation numbers look as follows 



where summation indices n, m = 0, . . . , oo in both (|2.3I1 and (|2.4p formulae. 

In the X-operator representation, the single-site part of Hamiltonian (|2.ip can be written as 

^o = EEE^-^r~' (2.5) 



where 



A„m = nie — /i) m(e' — \i) + + TO)(n + m — 1). (2.6) 



Terms describing an inter-site transfer in Hamiltonian (j2.ip are transformed in a similar way. 

Our primary goal is to calculate the two-time temperature boson Green's functions and 
((clc"*")), which describe an excitation spectrum and make it possible to investigate the conditions 
of the system's instability with respect to the spontaneous symmetry breaking and the appearance 
of a BE condensate. As follows from definitions 



nm rs 

((cdO). - V^^VITT((xr'™^-™+i|x;-^+i^-))^ . (2.7) 



nm rs 



We will use the equation-of-motion method for the evaluation of X-operator Green's functions. 
For the first one, from relations (|2.7p one could write 

+ . (2.8) 
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Let us write the commutators 

T ^n,m:n+l,m £7- -i / ■> •> \ ^n,m;n+l,m /r, (\\ 

[X, ' ' ^ , bll = (5/^ Vn- + 1 (^X; ' - X; ' ^ ' j , (2.10a) 

r -vn,m;n+l.m 71 r / / I FT -\An,m;n4-2.m / — x^n — l,m:n+l,Tn \ /o inu\ 

[X;' ,bi\^ 6iiiVn + 2Xi' -V"'-'^; ), (2.10b) 

[X, ' ' ^ ] = Su WmXi ^ • - ^/m + lXi ^ - ^ J ^ (2.10c) 

V^i ,CjJ = Oh ( + IXj ' -^/mX^' '1. (2.10d) 

The latter are originated from the commutation of an initial X-operator with the inter-site transfer 
terms of the Hamiltonian, thus producing the higher-order Green's functions 

{{xrb,\x;+^'^-'^^^)U , {{xrb+\x;+^''^-'^^^))^ (2.11) 

where Xf' stands for operators on the right-hand side of expressions (j2.10ap - (j2.10d[) . 

Decoupling of functions (12.111) in the random phase approximation (RPA) is performed in the 
following way: 

{{xrb,\x;+'^^--^n)u « {xnmx;+'^^--^n)^ + {b,){{xr\x;+'^^--^n)u. . (2.12) 

In the case of the normal phase (which will be studied herein) (bj) — (6^) — 0. Thus, retaining 
only the averages {Xi") of diagonal X-operators we have 

j^„,™;„+l,™^ ^] « A„„,Xr'™^"+l^" + V^Qnm ^ ^ub, + V^Q,„n ^ ^IjCj (2.13) 

3 3 



and equation (|2.8I) can be rewritten as 



^ r r r Qr. 



27r fioj — Ar, 



^ty((6,|x;+i.--)). + ^g±I^^t;;.((c,|x;+^— )). . (2.14) 



J 3 



The following notations are introduced 

^ ; Ti^n,m:n,m _ ^n+l,in\n+l,m\ * \ (O 

for the occupation difference of adjacent levels and the related transition energies when the number 
of Bose particles in the ground state (with the energy e) on the site increases by one. 

Proceeding from X-operators in equation (|2.14p to the Bose operators b and 6+ according to 
definition (j2.3p we obtain 

((^/I^P»- = ^^^p5o(u;)+5o(^)(E*y((^.l^p))- + E^y(fel^p))-)' (2-16) 

j 3 

where the function 

5oH = Er%^(- + i) (2-17) 

nm 

has the meaning of the unperturbed Green's function for bosons residing in the single-site ground 
state. 

Equations of motion for "mixed" Green's functions ((c|6+)) are obtained in the way similar to 
the above described scheme. Using decoupling p.l2p one can write 
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which resuhs in the equation 

((Q|6;)).=5oM(E*!i((^^-|^p))-+E*^.((c^l^p))-)- (2.19) 
Here, similarly to and (PT7)) 

r^l / -^n,m;n,m -^^n,m+l;ri,m+l\ a/ \ •> I I \ \ Qnm / , -i \ 

^7un — - /; ^nm — ■^n,7n+l ~ ^n,m , ffol^j — / , ^ ^, [^^+^), 

nm 

(2.20) 

and the function (7o(^) unperturbed Green's function for bosons residing in the excited state. 

By means of the Fourier transform 



N 
1 



one can proceed to the momentum representation obtaining a system of equations 

h 



^/n I Ij; I -i- (lr\ \ (jJ { ) ) . . H- (lr\ \ 111 I //' ^ 

27r- 

= 5o(^)i;'((^l&+)),,^ + g',W,{{c\h+)),.^^ , (2.22) 



where tq , t'^ and t'^ stand for the Fourier transforms of hopping parameters. 
A pair of equations for other Green's functions are obtained in a similar way 

((c|c+)),,. = Ag^(^) + g^(c.)t;'((6|c+)),,^ + g',{Lo)t'^{{c\c+)),^^ . (2.23) 
Solutions of equations (|2.22p and (|2.23l) are as follows: 



((c|c+)),.^ = — -^-^5^(u;)(l - 50(^)^9), 

= :^-^5o(^)5o(^)^;' = , (2.24) 



where 



Dq{iu) = 1 ~ g^{iu)tq - g'^{iu)t'q + g^{io)g'^{io) [t^t'^ - (t^')'] • (2.25) 

The equation Dq{Ld) = gives the excitation spectrum which is obtained here in the RPA. On 
the other hand, the divergence of boson Green's functions (I2.24p at the zero values of wave vector 
and frequency is the criterion of instability with respect to BE condensation [s^ , thus giving 
the following condition 

i?,=o(c^ = 0) = 0, (2.26) 
which can be rewritten in the explicit form 

1 - ga[u)tq - g',[u)t'^ + gaW,[u) [tqt'^ - (t^')'] - 0, (2.27) 

where 

5o(0) = -Er-%7 (" + !)' ^ (™ + l)' (2-28) 

^-^ (n + m)U — Li ^-^ (n + m)u +0 — 11 

nm nm 

and 5 = e' — e is the excitation energy. 
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We should point out that divergence of the ((fc|t'^))o.o function correlates with the appearance 
of the BE condensate in the ground state while at the divergence of the ((c|c+))o,o function, BE 
condensation takes place in the excited state. In general, both condensates appear simultaneously 
except the case t"^ = (e.g. due to symmetry reasons) where these effects become independent and 
only the one type of condensate arises in the instability point. 

Equation (j2.27p . mutually relating the chemical potential, hopping parameters and temperature, 
allows us to construct spinodal surfaces (or lines) in the above mentioned coordinates and to find 
the temperature of the phase transition to the phase with BE condensate (so-called SF phase) 
where such a transition is of the second order. Below, this problem (especially the issue of the 
phase transition order) will be investigated more in detail. 



3. NO phase instability in HCB limit 

Let us consider now a simple special case of the HCB limit when occupation numbers in the 
|ri, m) state are restricted by a condition n + to ^ 1. In the framework of the model, it formally 
means U ^ oo. 

In this case, the model becomes a three-level one with the local energies 



Aoo — 0, 

and the following transition energies 



An 



Ai 



A' 



00 



Thus, equation (|2.27|) can be rewritten in the form 



Q 



00 



00 



where 



Qoo — 



to - 



1 - e'^^ 



QooQ'o 



00 



[tot'o - {t'oY 



0, 



Qoo ^ 



1 



I _(_ Q^fJ. _|_ g/3(Ai-5) 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



in the zero approximation with respect to hopping. 
We take into account (according to esti- 



mations made in |l5l . |25|) that boson wave 



functions in adjacent potential wells overlap in 
greater extent in the excited states compared to 
the ground ones. Accordingly, we shall put here 
to — 0. For a centrosymmetric lattice and in 
the case of different parity of wave functions of 
ground and excited states we have also = 0. 
Finally, we follow a usual convention of the BH 
model for optical lattices taking < 0. In this 
way equation p.3p can be reduced to 



® 
re 

CL 

E 

CD 
1- 



\t'n\ 



1 



5 - ^1 + ef^t' + c'3(a'- 



1. 



(3.5) 




Its solutions determine the stability region 
boundaries of the normal (NO) phase. Respec- 
tive lines of spinodals are numerically calcu- 
lated and presented in figure [T] (here and below 
the energy quantities are given in units of |tg|). 

As illustrated in figure [TJ at (5 < jipl 
spinodals surround an asymmetric area in the 



-0.6 -0.3 0.0 
Chemical potential ^ 



Figure 1. Lines of the NO phase instability 
(spinodals) with respect to the appearance of 
BE condensate in the {Q, fi) plane in the HCB 
limit at various excitation energies {to = 0, 

\t'o\ = 1, 4 = 0). 
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(8,/i) plane which is located between the points /i = 5 — I^qI and ^ = of the abscissa axis. In 
this region, the NO phase is unstable; this is connected with the appearance of BE condensate. At 
S < |tQ|/2 and /i > the backward path of spinodal is observed and a lower temperature of the 
NO phase instability appears, thus suggesting a possibility of the SF phase existence in the inter- 
mediate temperature range (so-called "re-entrant transition"). However, as will be shown further, 
in the mentioned region a real thermodynamic behaviour is even more complicated. The order of 
the NO-SF transition can change to the first one and the SF-phase remains stable up to the zero 
temperature. 



4. Phase diagrams in MFA 

For a more detailed treatment of the NO-SF transition issue, let us study the thermodynamics 
of the considered system in the HCB limit, thus reducing the problem to a three-state model with 
the Hamiltonian 

h=y: \xr + E t.^xi'^x^' + E + E tu^i'^r + x^'^n (4-1) 

ip ij ij ij 

where the shorthand notations are used 

|0) = |00), |1) = |10), |2) = |01); Ao = Aoo, Ai = Aio, A2 = Aoi . (4.2) 

Possibility of BE condensation will be studied in the MFA. Average values of creation (annihi- 
lation) operators for Bose particles in the ground or excited local state 

= = (Xf ) (^ (6.)), e = {XD = {XD (^ (q)) (4.3) 

play the role of order parameters for the SF-phase. Hence, the mean-field Hamiltonian is as follows: 

ip 

+ E [tovixi^ + xr) + uixr + xf) + mxr + xf^) + f^vixf' + x^] . (4.4) 

i 

Self-consistency equations for parameters 77 and ^ 

7? = Z-iSp[XiOexp(-/37?MF)], e = ^"'Sp[Xf exp(-/37?MF)] (4.5) 

are equivalent to the condition of minimum of the grand canonical potential Q — —QlnZ, where 
Z = Spexp(-/3i7MF)- 

Limiting our consideration to the case of particle hopping only through excited states (ip 7^ 0, 
to — tQ — 0) we can diagonalize Hamiltonian (|4.4p by a rotation transformation 



(4.6) 



where 

cos2^= , sin2^= , ^'^"'^ (4.7) 

v/(A2-Ao)2+4(t[,0'' v/(^2-Ao)2+4(t[,e)2 

and A2 — Aq = (5 — /i. In terms of operators X'''^ = \ f){s\ 

HMF = N\Qe+J2~XpXfP. (4.8) 




cos?? 





— sin d \ 


( |0) 





1 






sin I? 





cos I? J 


V |2) 
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New energies of single-site states are 



(5 — /i 




Ai 



In the new basis 



-(Xf - XOO) sin 21? + (Xf - cos 21?, 



which yields after averaging 



(4.9) 



(4.10) 



(4.11) 



Taking into account that (X^p) — Z ^ exp(— /3Ap), Z ~ exp(— /3Ap) we come to the equation for 
the order parameter ^: 



-/3Ao 



^(J_^)2+ 4(^/^)2 V 

(4.12) 

Solution ^ = corresponds to the NO phase. A 
nonzero solution describing the BE condensate 
is obtained from the equation 



1 



= 1. 



Z ^(5_^)2+4(t^^)2 V 

(4.13) 

In the limit ^ — ?> this equation determines the 
line where the order parameter for the SF phase 
tends to zero. One can readily see that it coin- 
cides with spinodal equation p.Sp thus defining 
the line of the second order NO-SF phase tran- 
sition (when just the transition of such an order 
takes place). 

Numerical solutions of equation (|4.13p make 
it possible to study the behavior of the order 
parameter ^ depending on chemical potential \x at 
the main, at negative values of chemical potential 
transition to the SF phase is of the second order. 




-0.5 0.0 
Chemical potential ^ 

Figure 2. Dependences of the order parameter 
^ on the chemical potential /i for the reduced 
three-level (HCB) model at various tempera- 
tures indicating the possibihty of the first or- 
der phase transition at low enough temperatures 

(5 = 0.1, \t[A = \). 

various temperatures as illustrated in figure[5] In 
the parameter ^ changes smoothly and the phase 
But at /J, > and low enough temperatures, the 




Chemical potential n Chemical potential n 



Figure 3. Low-temperature behaviour of the order parameter ^ for the reduced three-level (HCB) 
model at zero and negative excitation energies 5 and various temperatures (jfo| = 1). 
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Figure 4. Lines of the NO-SF phase transi- 
tion in the (0, jj.) plane at various excitation 
energies S (|fo| — 1). 



Figure 5. An illustration of discrepancy between 
the spinodal curve and the real line of the first 
order phase transition for 6 = 0.01 (|to| = 1). 



^{^J-) dependence has an S-like bend. In this case, the first order phase transition with an abrupt 
change of the parameter ^ takes place. This phase transition occurs at a certain value of the 
chemical potential which could be calculated using the Maxwell rule or considering the minimum 
of the grand canonical potential 51 (/i) as a function of the chemical potential (see below). Obviously, 
the point of ^ nullification does not anymore correspond here to the phase transition. 

Similar behaviour of the parameter ^ holds even at zero excitation energy (i5 = 0) where the 
first order phase transition remains for nonzero temperatures whereas at T = its order changes 
to the second one (figure [3| ■ At negative values of S (which corresponds to inversion of e and e' 
levels and to hopping between ground states) the second order of the transition is preserved in the 
low-temperature region close to T = transforming to the first order transition at the temperature 
increase and recovering henceforth (figure [3]) . 

Changes in the NO-SF phase transition order and localization of the corresponding tricritical 
points are depicted in figure 21 where phase diagrams are given for various values of the excitation 
energy S. At temperatures lower than tricritical, spinodal lines and phase transition curves come 
apart as one can see comparing figures [T] and 01 At small values of S, the discrepancy is quite 




Figure 6. Appearance of two tricritical points 
at zero and negative values of excitation en- 
ergy S {\t'o\ = 1). 



Figure 7. Lines of the NO-SF phase transition 
in the (|to|,/i) plane at various temperatures O 
(energy quantities are given in units of 5). 
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significaiit (figure [5]). In the case of (5 < 0, two critical points appear at a certain distance; tfie 
latter tends to zero at 6 = ScHt ~ ~0.12|iQ| and the first order phase transitions at a further 
increase of S (figure IH]) is suppressed. 

Phase diagrams in the (|tg|,//) plane at various temperatures for (5 > are depicted in figure [7] 
with indication of tricritical points. In distinction to the standard two-level HCB model [34] (where 
the SF phase transition is of the second order) the diagrams are asymmetric. In the limit T — > 
for n > the first order transition occurs at /i = (v^ — y^ji^)^ (see the next section) whereas for 
/X < they are of the second order on the line n = 6 — IIq]. 



5. Phase separation at fixed boson concentration 

Let us consider now the thermodynamics of the model at a fixed concentration of Bose particles. 
We will utilize a connection between the concentration and the chemical potential of bosons which 
can be established using its definition in such a form 

n^{n'l + n^)^{Xl'+Xr) (5.1) 



or basing on the relationship 



d{n/N) 



(5.2) 



In the first case similarly to equality (I4.10p one can obtain a relation 

Xl^ + Xf = X," + sin^ d + If cos^ d + (Xf + Xf ) sin ^ cos ^ (5.3) 
which results in 

n = {Xl^) + (X™) sin2 d + (Xf ) cos2 d = 



1 



6 — jJL 



2 2^{5-iiY+^{UY 
In the second case, taking into account that 



1 



S — fi 



2 2V(<5-/.)2 + 4(i;,0^ 



-l3\2 



Z ^ c^^ + e" 



5-/1 ^ 
2 / 



(5.4) 



(5.5) 



and differentiating with respect to /i, one can come to the same expression as (j5.4p . 

There are different relationships between n and in NO and SF phases; in the last case, a 
nonzero value of ^ (a solution of equation (j4.13p 'l should be substituted into expression (|5.4p . Order 
parameter ^ has a jump at the first order phase transition, so a stepwise change of concentration 
n takes place. In the n = const regime (at the value of n in the region of step) it means a 
phase separation into two phases with different concentrations: the NO phase (i^ = and a larger 
concentration of bosons) and the SF phase (^ ^ and their smaller concentration). 

The above described situation is illustrated in figure[51 where the numerically calculated (O, n) 
phase diagrams are presented. At (5 > 0, phase separation region spans up to tricritical tempera- 
tures. When S goes to zero and finally reverses its sign, the shape of the separation region changes 
in a peculiar way moving off abscissa axis (figure [S]) . Now the phase separation begins at nonzero 
temperatures and vanishes at 5 < Scrit- the line of the second order phase transition remains 
only. At the further increase of \S\ (in the /i < region) the (0,n) diagram becomes more and 
more symmetric, approaching by its shape the diagram known for the usual HCB model [3^ (see 
also |36|). 

Phase diagrams in the (|iQ|,n) coordinates are given in figure [H] where the regions of NO, SF 
and separated phases are shown at various temperatures. 
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Figure 8. Lines of the NO-SF phase transition and the phase separation region in the (0,n) 
plane at various excitation energies 5 including the case of small, zero and negative values of S 
{\t'o\ = l)- 



The case of the zero temperature can be studied more in detail in a pure analytic way. In this 
limit there are three branches of order parameter ^ as a function of the chemical potential (see 
figure [2): 



(3):^-0. 

After elimination of ^ parameter, one can obtain the grand canonical potential fl as follows: 

(2): ^l/N=iS/\t',\-l)^^, 
'O, /i<0. 



(5.6) 



(3): n/N : 



-/X, jj, > 0. 



(5.7) 
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Figure 9. Phase diagram with the indication of possible phases (above) and lines of the NO-SF 
phase transition in the (|toli A*) plane at various temperatures (energy quantities are given in 
units of S). 
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Differentiating expressions (j5.7p with respect to /i we have 

1 fi-S 



(l):n 



2 2\t',\ 



(3):n=^' ^-"^^ (5.8) 
\l, /i>0. 

At the first order phase transition from the SF phase to the NO phase, the order parameter ^ 
jumps from branch (1) to branch (3). This occm's at the /i = /i* = i\/\t[)\ — VS)^ value given by 
equality of respective grand canonical potentials ^(i) — ^{3)- Then boson system separates into 
SF and NO phases with concentrations of bosons: 

1 fi* -S 

UNO = 1- (5.9) 



6. Discussion and conclusions 

As was shown in this work, the transition to the SF phase (the phase with BE condensate) in 
the Bose-Hubbard model with two local states (the ground and excited ones) on the lattice site 
can be of the first order in the case, when the particle hopping takes place only in the excited 
band. Calculations and estimates for optical lattices give evidence of significant distinction be- 
tween hopping parameters to and ti in the ground and excited bands, respectively. It follows from 
estimates [25] that ii/tg ss 30 — 50 depending on depth Vq of local potential wells (one can produce 
effect on Vb changing the intensity of laser beams which create an optical lattice). Similar results 
are obtained in the studies of quantum delocalization of the adsorbed hydrogen atoms. One can 
see from calculations 22 , 2^ of energy spectrum of the H-atom subsystem on the nickel surface 



that the ground-state band has a negligible bandwidth. At the same time, for excited bands, the 
bandwidth varies in the range from 15 to 45 meV (depending on the excited state symmetry and 
on the crystallographic orientation of metal surface) , being mostly of the order of half the corre- 
sponding excitation energy Asa = ^ £o ■ There are, however, the cases of strong delocalization 
(e.g. H on the Ni(llO) surface) where the excited bands overlap, and the width of the lowest one 
is of the same order as Asa [22j] . 

The values of hopping parameters greatly increase at the decrease of Vq ; the distance between 
the local energy levels becomes smaller in this case (see (stI. [38|). It is one of the possible ways 
of changing the relation between the hoppingparameters and excitation energy (|to| and 5 in 
our model). Another possibility (discussed in [39|) is connected with an essential reduction of the 
energy gap between local s- and p-levels due to sufficiently strong interspecies Feshbach resonance 
in the presence of Fermi atoms added to the Bose system in optical lattice. 

Along with investigation of BE condensation in the excited band (or bands PxiPy iPx,Py,Pz) 
in two- (three-) dimensional case) on condition that certain concentration of Bose-atoms has been 
created in the band by optical pumping [25, 38], an attempt was made in (40| to study the effect of 
excited bands on the physics of BE condensation in the lowest (s-) band (when the s-band hopping 
is taken into account). The case of finite values of the one-site interaction U was considered. 
The possibility of the re-entrant behaviour of the MI-SF transition was claimed. However, the 
order of phase transition was not investigated; the consideration was restricted to the case of zero 
temperature. As we show in this work, re-entrant type dependence on T or /Lt takes place only for 
spinodals and the return to the initial MI phase from the SF phase could be possible only in the 
case of the second order phase transitions. In reality, the order of phase transition changes to the 
first order in this region. In the HCB limit (no more than one particle per lattice site), it takes 
place mainly at positive values of chemical potential of particles; at /i < 0, the transition remains, 
for the most part, of the second order. The region of existence of SF phase is restricted, as a 
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whole, to the interval — |to| < fi < Itgl, while excitation energy should obey the inequality S < \tQ\. 
We have constructed the corresponding phase diagrams and established localization of tricritical 
points, where the order of phase transition changes. The separation on SF and NO phases at the 
fixed particle concentration is investigated; the conditions of the appearance of phase-separated 
state are analyzed. 

It should be mentioned that phase diagrams in figures [2Hni are close by their shape to the 
diagrams obtained in the framework of Bose-Hubbard model for Bose atoms with spin S* = 1 in 
optical lattices [2§|. The excited levels are formed in that case by the higher spin single-site states 
and corresponding interactions of the "ferromagnetic" or "antiferromagnetic" type (the Hund-rule- 
like splitting), while the hopping parameter is taken the same for all bands. The similarity of the 
mentioned diagrams points out to the fact that the role of the excited states in the change of the 
phase transition order in going to the phase with the BE condensate is the same in both cases. 
Distinction, however, consists in another genesis of the single-site spectrum. In our model, in the 
limiting case of HCB there are no effects connected with the level splitting due to interaction; the 
excited single-particle states are taken by us into account instead. 

The consideration developed in this work can be extended to the systems with the close or 
degenerate excited local levels. Generalization of the model by adding inter-site interactions is 
also important. It could even make it possible to take into consideration other phases (density- 
modulated or supersolid) besides NO and SF ones. 

We finally emphasize that the hopping parameter t'^^ in the excited band can be positive; in 
particular, this concerns the p-bands [39|. In such a situation, the condensation takes place into 
states with wave vector Q on the boundary of the Brillouin zone, while the order parameters (cq), 
(cq) describe the modulated condensate. Since t'g = — tg , the results obtained in this work are 
also valid (with jigl in place of Itgl) in that case. 
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ABOcraHOBa MOAeib 5o3e-Xa66apAa 

B rpaHMUi >KOpCTKI1X 6030HiB 

I.B. CracKDK, O.B. Be;iMHKO 

iHCTMTyT c()i3MKM KOHAeHCOBaHMx CMCTeM HAH VKpaiHn, By/i. I. CBeHuiubKoro, 1, 79011 /IbBiB, VKpaiHa 

floc;iiA>KeHO nepexifl y cfiasy 3 6o3e-KOHAeHcaT0M y flB030HHiM MOAe/ii Bo3e-Xa66apAa 3 nepeHOCOM 
HacTMHOK jiMLiie y 36yA>KeHiM 30Hi. Po3rjiflHyTO HecTiMKicTb, noB'H3aHy 3 TaKMM nepexoAOM, skum BiA6yBa- 
exbCH npn eHepriqx 36yA>KeHHH 5 < \t'g\ (p,e \t'g\ - napaMOxp nepecKOKy nacTMHOK). y rpaHnui >kopctkmx 
6o30HiB BnqBJieHD 3BopoTHiM xiA cniHOAa/ieM B o6;iacTi AOAaTHMx 3HaHeHb xiMiHHoro noTeHuia/iy. Bcxa- 
HOB/ieHo, mo piA 4)a30Boro nepexoAy y UbOMy BunaAKy 3MiHioeTbCfi Ha nepinuM i noBepHeHHn ao Hop- 
Ma/ibHOi c()a3M HacnpaBAi He BiASyBaeTbca. flepexoAn nepiiioro poAy ichymb i npn BiA'eMHiw eneprii 5 (3a 
yMOBM 5 > (Scrit ~ — 0.12|tg|). ripn ^ < 4)a30BMM nepexiA 3ajiMLijaeTbCH nepoBaxHo nepexoAOM APyro- 
ro poAy. npoaHa;ii30BaHoxiA napaMBxpa nopfiAKy 6o3e-KOHAeHcaTy, no6yAOBaHo (paaoBi (6,/^) i (liol,^) 
AiarpaMM, BCTaHOB/ieho ;ioKajii3auiio tpmkpmtmhhmx tohok. SHawAeHo yMOBw, npn hkmx a^i^ cfiiKCOBaHoi 
KOHueHTpauii nacTMHOK Hacxynae po3LijapyBaHHfi Ha Hopivia/ibHy (pasy i 4)a3y 3 6o3e-KOHAeHcaTOM. 

KnioHOBi cJioBa: Monenb 5o3e-Xa66apAa, xopcTKi 6o30hi/i, 6o3e-KOHfl,eHcau,isi, 36yM>KeHa 30Ha 
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